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ABSTRACT 

This paper describes recent work in the Raman lidar liquid water cloud measurement technique. The 
range-resolved spectral measurements at the National Aeronautics and Space Administration Goddard 
Space Flight Center indicate that the Raman backscattering spectra measured in and below low clouds 
agree well with theoretical spectra for vapor and liquid water. The calibration coefficients of the liquid 
water measurement for the Raman lidar at the Atmospheric Radiation Measurement Program Southern 
Great Plains site of the U.S. Department of Energy were determined by comparison with the liquid water 
path (LWP) obtained with Atmospheric Emitted Radiance Interferometer (AERI) and the liquid water 
content (LWC) obtained with the millimeter wavelength cloud radar and water vapor radiometer (MMCR- 
WVR) together. These comparisons were used to estimate the Raman liquid water cross-sectional value. 
The results indicate a bias consistent with an effective liquid water Raman cross-sectional value that is 
28%-46% lower than published, which may be explained by the fact that the difference in the detectors' 
sensitivity has not been accounted for. The LWP of a thin altostratus cloud showed good qualitative 
agreement between lidar retrievals and AERI. However, the overall ensemble of comparisons of LWP 
showed considerable scatter, possibly because of the different fields of view of the instruments, the 350-m 
distance between the instruments, and the horizontal inhomogeneity of the clouds. The LWC profiles for 
a thick stratus cloud showed agreement between lidar retrievals and MMCR-WVR between the cloud base 
and 150 m above that where the optical depth was less than 3. Areas requiring further research in this 
technique are discussed. 


1. Introduction 

Liquid water clouds play an important role in the 
earth’s energy balance by scattering and absorbing solar 
and terrestrial radiation (e.g.. Turner et al. 2007b). The 
influence of the clouds on the radiation balance depends 
on their microphysical properties such as liquid water 
content, number concentration and droplet size distribu- 
tion. These properties are affected by the microphysical 
properties of the aerosols from which they form, and the 
thermodynamic and dynamic environmental conditions 
(e.g., humidity structure, vertical velocity). For the study 
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of the effect of aerosols on clouds, it is important to 
measure the aerosol and cloud properties in the lower 
portion of the cloud; this being the region of droplet 
formation and rapid growth (e.g., Mordy 1959; Lamb 
and Verlinde 2011). Also, it is desirable to measure 
these properties in the same atmospheric volume using 
a single instrument because different instruments have 
different spatial and temporal resolutions and volumes, 
resulting in uncertainties in studying aerosol-cloud pro- 
cesses. The Raman lidar offers the potential to measure, 
simultaneously and in the same atmospheric volume, 
aerosols, water vapor, and liquid cloud properties in this 
critical droplet formation region. For that reason, we 
investigated further this measurement technique. This 
paper focuses on the liquid water measurement tech- 
nique since it has not been well studied and verified. 
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The first report of the range-resolved observation of 
the Raman scattering by liquid water cloud (fog) was 
Bukin et al. (1983). Since that time, Melfi et al. (1997) 
have discussed the use of Raman scattering from cloud 
droplets to quantify the liquid water content of clouds. 
They measured excess (apparently supersaturating) sig- 
nal in the Raman water vapor channel from the cloud and 
attributed it to the Raman scattering by liquid droplets. 
Whiteman and Melfi (1999) extended this work and 
presented a technique for computing the liquid water 
content, mean radius, and number concentration of 
droplets by using the Raman and Mie backscattering 
signals from clouds. However, the measurements were 
with broadband filters, where the liquid and vapor sig- 
nals were together in a single channel. Veselovskii et al. 
(2000) measured the Raman water vapor and liquid 
water backscattering separately, but not simultaneously, 
by using interference filters. They observed increases in 
the liquid water signals both in the clouds and in the 
boundary layer compared to that above the top of the 
boundary layer. They attributed the liquid water signal 
in the boundary layer to water-coated aerosols. Rizi 
et al. (2004) measured water vapor and liquid water 
signals simultaneously in separate channels and con- 
cluded that estimating the liquid water content is pos- 
sible when sampling low-level clouds if the spectral 
efficiencies of the lidar channel can be measured and/or 
estimated, the water vapor calibration constant is eval- 
uated, and the ratio between the differential Raman 
backscattering cross sections of water vapor and liquid/ 
droplet is known. Whiteman et al. (2010) demonstrated 
airborne retrievals of cloud liquid water content, droplet 
radius, and number concentration where the droplet size 
distribution obtained from those measurements were in 
good agreement with the literature. 

One major problem in the previous studies is that 
there were few comparisons of the lidar-derived cloud 
properties with those obtained with other measurements 
in order to attempt to validate the technique. Because 
the magnitude of the Raman backscattering cross sec- 
tion of water droplets is not well known (Veselovskii 
et al. 2002) and the wavelength of fluorescence of air- 
borne particles likely overlaps the Raman wavelength 
(Manninen et al. 2009), we first wanted to confirm with 
spectral measurements that the lidar does detect the 
Raman liquid water backscattering from clouds and that 
those measurements are in reasonable agreement with 
expectation. We then calibrated the lidar signal in an 
appropriate manner. 

Thus, in this study we first measured the backscat- 
tering spectrum from a water cloud over the National 
Aeronautics and Space Administration (NASA) Goddard 
Space Flight Center (GSFC) in Greenbelt, Maryland. We 
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Fig. 1. Differential backscattering cross sections of Raman LW, 
WV, and N 2 excited by 354.7-nm radiation at 300 K. Curve for 
LW uses an arbitrary scaling. Transmission curves (blue) of the 
IFs for the Raman lidar at the ARM site are referenced to the 
right axis. 


then analyzed data obtained with the Raman lidar at the 
Southern Great Plains (SGP) site of the U.S. Department 
of Energy’s Atmospheric Radiation Measurement Pro- 
gram (ARM), where calibration coefficients for the lidar 
were obtained by comparisons with 1) the Atmospheric 
Emitted Radiance Interferometer (AERI) and 2) the 
millimeter wavelength cloud radar and water vapor ra- 
diometer (MMCR-WVR) measurements. The content 
of this paper is as follows. Section 2 describes the methods 
for measuring liquid water with the Raman lidar. Sec- 
tion 3 shows the results of the spectral measurements 
acquired at NASA GSFC and time series of measure- 
ments from the ARM SGP site. Section 4 discusses the 
current uncertainties of this technique by focusing on 
aerosol fluorescence and the Raman cross section. 
Section 5 summarizes this study. 


2. Raman liquid water measurement technique and 
instrumentation 

a. Principle of Raman liquid water measurement 
technique 

The measurement of liquid water content by lidar is 
based on Raman backscattering by cloud droplets. The 
laser emits pulses of light (we use the tripled Nd:YAG 
wavelength of 354.7 nm) and light backscattered by 
cloud droplets is collected with the receiving telescope. 
The lidar return contains the different backscatter 
photons that are separated into Raman (403.2 nm for 
liquid water, 407.5 nm for water vapor, and 386.7 nm for 
nitrogen) and Rayleigh-Mie (354.7 nm) components 
using dichroic beam splitters and interference filters 
(IFs). Figure 1 shows the differential Raman backscat- 
tering cross sections of these atmospheric components 
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Table 1. Specification of lidar spectrometer at NASA GSFC. 



Transmitter 

Laser 

Tripled Nd:YAG 

Wavelength 

354.7 nm 

Pulse energy 

260 mJ 

Pulse repetition frequency 

50 Hz 

Beam diameter 

3 cm (0.125 mrad divergence) 

Bandwidth 

~1 cm -1 


Receiver 

Primary telescope 


diameter and f-number 

41 cm, f/10 

Focal length 

0.75 m 

FOV 

0.5 mrad 

Rejection at 354.7 nm 

>10 6 

Spectrometer 

Czerny-Turner 

F /# 

f/9.7 

Slit width 

10 fj . m to 2.5 mm (adjustable) 

Grating groove density 

300, 1200, and 2400 line mnC 1 

Spectral range 

180-850 nm 

Resolution 

0.261,0.078, and 0.037 nm 

Detector 

ICCD camera 

Pixels 

1024 X 256 


and the transmission curves of the interference filters in 
use in the lidar at the ARM SGP site. 

b. Lidar spectrometer at NASA GSFC 

To confirm that the lidar detects the Raman liquid 
water backscattering and to study the possible inter- 
ference of fluorescence in the Raman liquid water mea- 
surement, we have implemented a lidar spectrometer 
at NASA GSFC for measuring the Raman spectrum 
backscattered by the atmosphere. Table 1 and Fig. 2 show 
the specifications and a schematic of the lidar spectrom- 
eter. It emits laser pulses of light at 354.7 nm with 
260 mJ pulse -1 with a repetition rate of 50 Hz. The light 
backscattered by the atmosphere is collected with a 
41-cm-diameter telescope (LX200EMC, Meade, United 
States). The light collected with the telescope is sepa- 
rated into the spectral component longer than —380 nm 
(transmission) and shorter than that (reflection). The 
transmitted light passes a blocking (shortcut) filter to 
suppress the strong backscattered light entering the 
fiber-optic bundle, which can induce fluorescence. The 
total rejection factor of the Rayleigh-Mie signal ex- 
ceeded 10 6 . The light is collimated using a lens and en- 
ters the fiber-optic bundle (Ceram Optec 19XUV200/ 
220P/2.0M) that possess a numerical aperture of 0.22. 
The bundle converts the shape of the light beam from 
circular to rectangular. The rectangular light beam en- 
ters the slit of a spectrometer after passing the lens to 
match the aperture ratio (f/value) of 9.7 of the spec- 
trometer. The spectrometer is a Czerny-Turner spec- 
trograph (Shamrock SR-750-A, Andor Technology, 


United States) equipped with a fast gated intensified 
charge-coupled device (ICCD) camera (iStar DH720, 
Andor Technology). The highest spectral resolution of 
the spectrometer is 0.261 nm when using a grating with 
a groove density of 300 lines mm -1 , as was the case for 
the measurements here. The actual resolution used for 
the measurements made here was — 1.6nm, realized by 
widening the entrance slit to increase the signal-to-noise 
(S/N) ratio of the measurement. The system measures 
the spectrum of the backscattered light at one height in- 
terval at a time through gating of the ICCD camera. We 
used a gate width of either 5 or 10 /xs, corresponding to 
height intervals of 750 or 1500 m, respectively. We accu- 
mulated the data for 5 min to increase the S/N ratio. In 
addition to measuring the spectrum, we measured si- 
multaneously the parallel and perpendicular components 
of the backscatter signal at 354.7 nm with respect to the 
polarization plane of the laser with photomultiplier tubes 
(PMTs) (R1924, Hamamatsu, Japan) and transient re- 
corder (Licel TR-20, Germany) to obtain the vertical 
distribution of the cloud. The depolarization ratio can be 
utilized to distinguish spherical particles from irregularly 
shaped ones. In this experiment, we used only the parallel 
channel because we did not optimize the alignment of the 
polarization of the receiving optics. The location of the 
lidar is Greenbelt, Maryland (38.99°N, 76.84°W and 
50 m above mean sea level). 

c. Lidar system at ARM SGP site 

The description of the Raman lidar system at the U.S. 
Department of Energy ARM site in the Southern Great 
Plains has been given elsewhere (Goldsmith et al. 1998; 
Turner and Goldsmith 1999, 2005; Newsom 2009), so we 
only briefly describe it here. Table 2 shows the specifica- 
tions of this lidar system. It emits laser pulses of light at 
354.7 nm and detects the backscattered light with a 61-cm- 
diameter telescope and 10 photomultiplier detectors. 
The liquid water channel was added to the system in Oc- 
tober 2005. For the detection of Raman liquid water 
backscattering, a broadband IF with the center wavelength 
(A C w) at 403 nm and full width at half maximum (FWHM) 
of 6.5 nm is used (Fig. 1). For the detection of Raman 
scattering from water vapor and nitrogen, narrowband 
IFs with A cw = 407.5 and 386.7 nm and FWHM = 0.3 nm 
are used (Goldsmith et al. 1998; Turner and Goldsmith 
1999). The ARM Raman lidar was upgraded in 2004 to 
use Licel detection electronics, where the vertical reso- 
lution of the data is 7.5 m vertically and 10 s temporally 
(Turner and Goldsmith 2005; Newsom et al. 2009). The 
acquired raw data are available from the ARM Data 
Archive (http://www.archive.arm.gov). The location of 
the lidar is near Lamont, Oklahoma (36.61°N, 97.49°W 
and 311m above mean sea level). 
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d. Method for computing liquid water mixing ratio 
from lidar data 

The standard approach to using the Raman lidar 
technique to derive liquid water content has been de- 
scribed previously (Whiteman and Melfi 1999; Rizi et al. 
2004). Assuming a constant enhancement due to the 
spherical configuration of the droplet (Veselovskii et al. 
2002), the ratio of the liquid water Raman signal and 


a molecular reference can be used to provide the liquid 
water mixing ratio in a direct manner (Whiteman and 
Melfi 1999). There are two potential contributions of 
signal in the liquid water channel that are not due to 
liquid water scattering: 1) vibrational-rotational (Raman 
scattered) lines of water vapor and 2) aerosol fluores- 
cence. Here, we introduce correction factors for them. 
The liquid water mixing ratio (LWR) at height z is then 
computed from the lidar data using 


LWR(z) =/. w /. 


^lw( z ) *av vAvv 


(z) 




LW-'WV 


( Z ) 


-// 


AE AZ, 


CB 


^Lw(z / ) ^WV^WV 


(z') 


Z„„-AZ„ 


PM 


dz', 

( 1 ) 


where P x (z) is the background subtracted signal in- 
tensity from the height z and the subscripts refer to the 
Raman liquid water (LW), water vapor (WV), and ni- 
trogen (N 2 ) channels; t x is the atmospheric (molecular 
and particle) transmission between the lidar and z at the 
wavelength x. We compute the molecular transmission 
from the radiosonde-derived atmospheric density and 


the particle transmission from the nitrogen Raman sig- 
nal and/or Rayleigh Mie signal using the method de- 
scribed by Sakai et al. (2003). We use the extinction 
coefficient derived from the Raman N 2 signal where the 
uncertainty in the derived extinction coefficient is lower 
than 100% or is estimated from the backscattering co- 
efficient () 3 P ) derived from Mie-Rayleigh and Raman N 2 
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Table 2. Specification of Raman lidar at ARM SGP site (after Newsom 2009). 


Transmitter 

Tripled Nd:YAG 
354.7 nm 
300-400 mJ 
30 Hz 

13 cm (—0.1 mrad divergence) 

—2 cm 1 

Receiver 

61 cm, f/9.3 
30%-40% 

2 mrad 
0.3 mrad 

Simultaneous photon counting and analog 7.5-m 
range resolution 


Wavelength 


Detection channels 

Center wavelength (CWL, nm) 

FWHM (nm) 

Unpolarized elastic (WFOV) (nm) 

354.7 

0.3 

Beam parallel elastic (NFOV) (nm) 

354.7 

0.3 

Beam perpendicular elastic (NFOV) (nm) 

354.7 

0.3 

Water vapor (WFOV and NFOV) (nm) 

407.5 

0.3 

N 2 (WFOV and NFOV) (nm) 

386.7 

0.3 

Temperature (NFOV) (nm) 

353.3 

0.25 

Temperature (NFOV) (nm) 

354.3 

0.25 

LW (NFOV) (nm) 

403.2 

6.5 


Laser 

Wavelength 
Pulse energy 

Pulse repetition frequency 
Beam diameter 
Bandwidth 

Primary telescope diameter and f-number 

Filter transmission 

FOV 

Wide (WFOV) 

Narrow (NFOV) 

Detection electronics 


signals by multiplying the extinction-to-backscatter ratio 
(we assumed 50 sr when the backscattering ratio is lower 
than 3 for aerosols and 20 sr when the backscattering ratio 
is equal to or larger than 3 for clouds). The coefficients f x 
are 1) the calibration coefficient of the liquid water mix- 
ing ratio (/ LW ), 2) the calibration coefficient of water 
vapor mixing ratio (/wv), 3) the fractional contribution of 
the Raman water vapor signal (Cwv) in the LW channel, 
and 4) the fractional contribution of aerosol florescence 
(/ae) to LWR; zcb is the cloud-base height and Az CB = 
300 m. The/ wv is obtained by a least squares fitting of the 
lidar data between 0.3 and 4.0 km in altitude to that ob- 
tained with the coincident radiosonde launched approx- 
imately 80 m from the lidar. The uncertainty in /wv is 
estimated to be less than 5% from the standard error of 
the least squares fitting of the lidar to the radiosonde data. 

The / LW includes the ratio of the effective Raman 
cross section of water vapor to liquid water and the 
detection efficiency of the receiving optics that is ex- 
pressed by (Whiteman 2003) 


da 


WV; 


/lw 




swv 


d\ 


da 


LW / LW dAXG 


( 2 ) 


da 


where da x /dfl is the differential backscattering Raman 
cross section at an exciting wavelength of 354.7 nm, f^is 
the detection efficiency, and G is the gain enhancement 


factor of the Raman liquid water cross section due to 
resonances of the spherical droplet (Veselovskii et al. 
2002). We assume that the beam overlap factors are 
equal for the two channels above 0.3 km in altitude 
based on the agreement of the lidar-derived water vapor 
mixing ratio with radiosonde without overlap correc- 
tion. For the purpose of this study and based on the 
range of temperatures (—30 K) encountered, the tem- 
perature dependence of the effective Raman cross sec- 
tion is negligible because the variation is less than 0.3% 
and 2% for the liquid water and water vapor, respec- 
tively. The aerosol fluorescence is estimated from the 
average of the uncalibrated LWR below the cloud base 
between z CB and z CB — Az CB , a region where no sig- 
nificant liquid water signal is expected in the absence of 
precipitation. Fluorescence is a spontaneous emission of 
a photon by an atom or molecule after a transition into 
an excited state due to absorption of the incident radi- 
ation. The decay time of the aerosol fluorescence is on 
the order of 10 ns (Pan et al. 2007), which is shorter than 
the resolving time of the lidar (50 ns). It should be noted 
that / LW , /wv - : and C W v are system constants, whereas 
/ AE can vary depending on the aerosol properties (size, 
shape, and chemical composition) in the cloud. We de- 
termined these coefficients for the lidar at ARM SGP by 
comparing with AERI and MMCR-WVR measure- 
ments as described in sections 3b and 3c. The theoretical 
values of / L w and C W v are calculated to be 0.081 and 
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0.035, respectively, from the Raman cross section of 
5.1 X 10 _3 ’m 2 sr _1 for liquid water (Ahmad and lies 
2001) and 6.9 X 10 _34 m 2 sr _1 for water vapor (Avila 
et al. 2004). The spectral shapes of the Raman band are 
given by Walrafen (1967) and Whiteman et al. (1999) 
for liquid water and Avila et al. (2004) for water vapor. 
The transmission ratio of the receiving optics of LW to 
WV is 1.2 from the manufacturer’s data on the beam 
splitters and interference filters used for separating and 
sampling these signals (Russo 2007). We assumed that 
G = 2 (Veselovskii et al. 2002). There are several pos- 
sible uncertainties in the coefficients as follows: 

1) The manufacturer’s sensitivity curves show that there 
is little difference between the quantum efficiency of 
PMTs at the liquid water and water vapor wavelengths, 
so these quantum efficiency values have been assumed 
to be the same. However, the sensitivity of individual 
PMTs can differ significantly and this possibility has 
not been accounted for. 

2) The Raman cross section of liquid water measured in 
laboratory experiments varies by 30% at an exciting 
wavelength of 354.7 nm (Ahmad and lies 2001; Li 
and Meyers 1990; Faris and Copeland 1997). 

3) Although likely of a much smaller magnitude than 
the two other effects mentioned, there still may be 
a question of whether the embedded dipole model 
simulations accurately account for all the spherical 
resonance of droplets (Veselovskii et al. 2002). 

The combination of these three could contribute a sig- 
nificant error to these theoretical calculations. It is for 
that reason that we will show later the calculation of 
these values using an optimization scheme. The liquid 
water content (LWC) is calculated as LWC = LWR X p, 
where p is the atmospheric density that was obtained 
from the radiosonde measurements. 

3. Results 

a. Backscatter spectrum measurement over NASA 
GSFC 

To study the spectrum of Raman liquid water signal 
from clouds, we measured spectrally resolved back- 
scattering from a low cloud over NASA GSFC from 
2013 to 2029 LT 15 November 2011. During the mea- 
surement period, a thin cloud was present at an altitude 
range between 1.8 and 2.0 km (top-right panels in Figs. 3 
and 4). Therefore, Figs. 3 and 4 show the spectra mea- 
sured in and below the cloud where the ICCD range gate 
was set to capture two different spectra: one between 
the altitude ranges from 1.6 to 2.4 km (Fig. 3), where the 
cloud was located: and the other between the ground 



Wavelength (nm) 


Fig. 3. Backscattering spectrum in cloud between 1.6 and 2.4 km 
in altitude over GSFC at 2024—2029 LT 15 Nov 2011. Thin gray 
curves show the theoretical Raman intensity fitted to the mea- 
sured spectrum with molecular number concentration ratio of LW 
to WV as 1 to 8.1. Top right inset shows the vertical distribution of 
backscattering intensity at 354.7 nm at 2026-2029 LT and hori- 
zontal dotted lines with arrow show the height range of spectrum 
measurement. 

and 1.6 km (Fig. 4), which was used as a clear-air spec- 
trum. The temperature in the height range of the cloud 
was between 278 and 280 K at 2000 LT based on the 
nearest radiosonde observation site (Dulles International 
Airport, which is a distance of 52 km away), indicating 
that the phase of cloud particles was likely liquid. 
Changes in the cloud structure during the measurement 
period were only a small concern because lidar mea- 
surements showed these variations to be small. 

The spectrum measured in the cloud (Fig. 3) shows 
that, although the signal-to-noise ratio of the data was 
low, the spectral intensity increased for the wavelength 
ranges from 400 to 420 nm and from 380 to 390 nm. The 
spectral shape in the 400-410-nm range fits well to the 
convolution of theoretical Raman cross sections of liquid 
water and water vapor with the molecular number 
concentration ratio of liquid water to water vapor (r) of 



Wavelength (nm) 


Fig. 4. As in Fig. 3, but below the cloud between 0 and 1.6 km in 
altitude over GSFC at 2013-2018 LT 15 Nov 2011 and with the thin 
gray curves fitted with a ratio of LW to WV as 1 to 31. 
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1/8.1. We determined this ratio by fitting the theoretical 
curve to the measured spectrum to minimize the func- 
tion defined by 

2 _ V ^obs( A i) “ ^Lw(^i) _ (1 ~ r )Aw(^i)] /o\ 

K or / \ %2 ’ WJ 

1=1 5/ obs( A i) 

where 7 obs , 7 LW , and 7 WV are the observed and theoret- 
ical Raman liquid water and water vapor intensities at A,- 
respectively. The uncertainty in the observed intensities 
(5 7 ibs was estimated by v/obs assuming the Poisson sta- 
tistics. We set Ai = 398 nm and \ N = 415 nm, with N = 
147, and exclude the data points between 406 and 408 nm 
because the observed signals might be saturated as dis- 
cussed below. 

The water vapor peak measured at 408 nm is 2.5 times 
lower than the theoretical value; this may be due to de- 
tector saturation. We estimate that the average LWC in 
the cloud (1.8-2. 0km in altitude) was 0.34 ± 0.08 gm" 3 
from the number concentration ratio of water vapor to 
liquid water in the spectrum (Fig. 3), Raman cross sec- 
tions of water vapor (Avila et al. 2004) and liquid water 
(Ahmad and lies 2001), and the radiosonde-derived 
vertically integrated water vapor content for the altitude 
range of the spectrum (1.6-2. 4 km in altitude), which is 
consistent with the value of low stratus cloud (Pruppacher 
and Klett 1997). The uncertainty in the LWC was esti- 
mated using standard error propagation from the un- 
certainties of 14% for the number concentration ratio of 
water vapor to liquid water; 10% for the Raman cross 
section of water vapor based on Penny and Lapp (1976); 
and 17% for the liquid water based on the standard de- 
viation of cross section reported by Ahmad and lies 
(2001), Li and Myers (1990) and Faris and Copeland 
(1997); and 2% in the radiosonde-derived water vapor 
content (Miloschevich et al. 2006). The large increase in 
the 380-390-nm range is the Raman backscattering by 
nitrogen (A cw = 386.7 nm). 

The spectrum below the cloud (Fig. 4) shows that the 
intensity increased in the spectral ranges from 400 to 420 
and from 380 to 390 nm. The spectral shape fits well the 
combination of Raman liquid water and water vapor 
spectra with r = 1/31 except for the peak at 408 nm, 
where the observed value was 3.2 times lower than 
theoretical value, possibly due to the detector satura- 
tion. It should be noted that the intensity at the 400-406- 
nm range (Raman liquid water region) is much lower 
than that in cloud (Fig. 3). The small increase in signal 
between 400 and 406 nm in Fig. 4 may be due to water- 
clad aerosols, as suggested by Veselovskii et al. (2000) 
and Rizi et al. (2004). From these results we confirmed 
that the lidar could detect the Raman liquid water 


backscattering from the clouds and distinguish the cloud 
from the cloud-free air in this case. However, it is de- 
sirable to improve the signal-to-noise ratio of the spec- 
trum to increase the confidence of the detection of the 
Raman liquid water signal. Kim et al. (2009) have also 
measured the Raman spectrum of the atmosphere at 
a wavelength range from 402 to 411 nm by using a mul- 
tichannel PMT. They found a difference in the shape of 
the spectrum between a cloudy day and a clear day, al- 
though there was ambiguity in the result because the 
spectral resolution of their measurement was too coarse 
to resolve the spectral shapes of the Raman bands. 

There is a possible enhancement at wavelengths lon- 
ger than 415 nm that could be due to fluorescence by 
aerosols in the air (we discuss this possibility later in 
section 4a). The width of the Raman nitrogen spectrum 
measured below the cloud (Fig. 4) is larger than that in 
cloud (Fig. 3). A possible explanation for that is the cross 
talk of the charge-coupled device (CCD) elements due 
to the strong signal from the lower altitude. 

b. Thin inhomogeneous clouds (altostratus) 
measurements over ARM SGP 

To obtain the lidar calibration coefficients (f x and 
C W v) shown in Eq. (1), we compared the liquid water 
path (LWP) obtained with the lidar and AERI over the 
SGP ARM site. The AERI measures the absolute in- 
frared sky spectral radiance between 3 and 19.2 pm di- 
rectly above the instrument with a resolution of 0.5 cm -1 
(Knuteson et al. 2004a,b). From the measured radiance, 
LWP and liquid cloud optical thickness were obtained 
using the mixed-phase cloud retrieval algorithm (MIX- 
CRA) developed by Turner (2005). The temporal res- 
olution of the AERI data was 19-25 s. The uncertainty 
(random error) in the AERI-derived LWP is less than 
4% for LWP < 50 g m~ 2 (Turner 2007). Importantly, the 
field of view (FOV) of the AERI is —46 mrad in FWHM, 
which is larger than that of the Raman lidar (0.3 mrad). 
The difference in FOV can cause errors in determining 
the calibration coefficients because of the sampling dif- 
ferences between the two instruments. For comparison, 
we also show LWP obtained with WVR. The WVR 
measures down welling radiant energy at 23.8 and 
31.4 GHz. The FOVs of the WVR are 5.9° and 4.5° (100 
and 79 mrad, respectively). From these observations, 
LWP and precipitable water vapor (PWV) are obtained 
using the microwave radiometer retrieval (MWRRET) 
algorithm developed by Turner et al. (2007a). The re- 
trieval approach uses an optimal estimation method- 
ology where the covariance of the observations is 
propagated through the retrieval to provide estimates 
of the uncertainty in the retrieved parameters. The 
temporal resolution of the WVR data was 20 s. The 
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Fig. 5. Temporal variation of LWP obtained with the Raman lidar (magenta), AERI (green), and WVR (blue) for the period 
0400-1000 UTC 9 Aug 2006 over SGP ARM site. Error bars show the standard deviation of data points averaged for 5 min for AERI and WVR, 
and the measurement uncertainty estimated from the uncertainty in the lidar signals using the Poisson statistics for the Raman lidar. 


uncertainty in the derived LWP is at least 20 gm - 
(Turner et al. 2007a). 

To compare the AERI and Raman lidar measure- 
ments, we needed to select clouds that had sufficiently 
low optical depth such that the lidar could profile 
through the entire cloud, since the AERI and WVR 
provide only LWP. Figure 5 shows the temporal variation 
of the LWP obtained for the period 0400-1000 UTC 
9 August 2006. The temporal resolution used was 5 min to 
improve the signal-to-noise ratio of the lidar retrievals. 
To compute LWP from the lidar data, we integrated 
LWC values vertically from below the cloud base to 
above the cloud top (3. 0-4. 5 km in altitude, see Fig. 7). 
The uncertainty in the LWP (LWC) is estimated from 
the uncertainty in the lidar signals assuming the Poisson 
statistics (Whiteman 2003). 

The qualitative agreement between the lidar and 
AERI is quite good, whereas the WVR showed values 
higher than these instruments for values lower than 
10 gm" 2 . We recall that the uncertainty of the WVR is 
approximately 20 gm -2 , which offers a possible expla- 
nation for the disagreement of the WVR results with the 
other instruments. We obtained calibration coefficients 


for the lidar [/ LW , C wv , and/ AE in Eq. (1)] by fitting the 
lidar data to the AERI-derived LWP by varying the 
values of /lw from 0.10 to 0.20 in intervals of 0.01, C W v 
from 0.01 to 0.03 in intervals of 0.01, and/ AE from 0 to 2 
in intervals of 0.1. We determined the values by per- 
forming a simple linear regression for each set of f x and 
Cwv values and finding the values that result in a slope 
that is closest to one, and an intercept that is closest to 
zero. Table 3 shows the results that provided the best fit 
to the data. The optimum values were determined to be 
/lw = 0.17, Cwv = 0-01, and/ AE = 1.0 for this case. The 
uncertainty in the f x values are estimated to be at least 
11% considering the statistical uncertainties of the lidar 
signals (<10%), and the measurement uncertainty in the 
AERI-derived LWP (<4%). However, the uncertainty 
should be larger than that because of the sampling un- 
certainty associated with the differing FOVs, the —350 m 
separation between the Raman lidar and AERI, and in- 
homogeneities in the cloud field. The fractional contri- 
bution of the Raman water vapor signal to the LW 
channel (i.e., the ratio of C wv Pwv to w) was 15% and 
the fractional contribution of the aerosol fluorescence to 
the LWR [i.e., the ratio of the second term to the first 


Table 3. Result of determining calibration coefficients of Raman lidar in this study. Asterisk denotes regression equations: LWP(Raman 
lidar) = slope X LWP(AERI) + intercept or LWC(Raman lidar) = slope X LWC(MMCR-WVR) + intercept. 


Instrument 

Date 

Parameter 




Regression* 


/lw 

Cwv 

/ 

Slope 

Intercept 

Correlation coefficient 

Number of data points 

AERI 

9 Aug 2006 

0.17 

0.01 

1.0 

0.97 

0.31 gm -2 

0.80 

41 

MMCR-WVR 

9 Apr 2008 

0.12 

0.03 

1.8 

1.00 

-9.7 X 10~ 3 gnT 3 

0.93 

12 

MMCR-WVR 

3 May 2008 

0.18 

0.02 

0.4 

1.00 

-6.1 X 10~ 3 gm~ 3 

0.96 

7 
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Fig. 6. Scatterplot of LWP obtained with AERI and Raman lidar 
for the period shown in Fig. 7. The equation in the figure is result of 
linear regression with y = LWP(Raman lidar) and x = LWP 
(AERI); r is correlation coefficient. Close-up for the small LWP 
values is the inset. 

term of the right side of Eq. (1)] is 41%, suggesting that 
the aerosol fluorescence was not negligible and that the 
derived LWR is sensitive to / A e for this case. 

A comparison of these values with the theoretical 
ones indicates that/ LW is 2. 1-2.2 times larger than the 
theory (/ LW = 0.081). We believe that the major reasons 
for the difference are the assumption that the detection 
efficiency of PMTs for LW and N 2 channels are equal 
and the uncertainty in the Raman cross section of liquid 
water as discussed in section 2d. We believe that un- 
certainty in G is smaller than the others because it is 
based on the electromagnetic theory and computation. 
The value of C W v is 3.5 times smaller than that expected 


by theory (C W v = 0.035). The possible reason for the 
smaller contribution of the Raman water vapor signal to 
the liquid water channel than expected for this case is 
that the cross talk of the liquid water channel differs 
from that computed from the characteristics of the op- 
tical components (e.g., spectral transmission of the di- 
chroic beam splitter for the WV channel placed before 
the LW channel). The value / AE = 0. 9-1.0 suggests that 
aerosol fluorescence intensity does not enhance in the 
cloud droplets and is almost the same inside the cloud as 
it is below cloud for this case. Using these coefficients, 
the lidar-derived LWP varied between 0 and 27 gm -2 
(Fig. 5). The correlation coefficient is 0.80 for 41 data 
points between 0400 and 0932 UTC (Fig. 6). The dif- 
ference in LWP between the lidar and AERI is large at 
0824 and 0906 UTC, possibly because of the inhomo- 
geneity of the cloud. Figure 7 shows the temporal and 
vertical cross section of backscattering ratio at 354.7 nm 
for the same time period as Fig. 5. The temporal reso- 
lution is 10 s and the vertical resolution is 60 m. We can 
see in the figure that the clouds with backscatter ratio 
higher than 10 were present between the altitude range 
of 3-4 km and the vertical distribution varied signifi- 
cantly with time. This variability likely was an influence 
on the differences in LWP between the lidar and the 
AERI. The lidar-derived LWP would be higher than that 
derived with AERI if the lidar pointed at the high LWP 
region of the inhomogeneous clouds, whereas it would be 
lower if the lidar pointed at the low LWP region. 

Figure 8 shows the temporal variation of the cloud 
optical thickness. The value for the lidar data was com- 
puted from the attenuation of the Raman nitrogen signal 
from an altitude of 3. 0-4.5 km — that is, the same altitude 
range as for the LWP retrieval. The lidar-derived value 
ranged from 0 to 2.43. It should be mentioned that the 
temporal variation of the optical thickness is similar to 
LWP, indicating the proportionality of the cloud optical 



Fig. 7. Temporal and vertical cross section of backscattering ratio at 354.7 nm for the same time period as in Fig. 6. 
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Fig. 8. As in Fig. 7, but for the optical thickness obtained with the Raman lidar (magenta) and AERI (green). Error 
bars show the standard deviation of the data points averaged for 5 min for AERI and the measurement uncertainty 
estimated from the lidar signals for the Raman lidar. 


thickness to LWP (Fig. 9). The regression analysis be- 
tween the lidar-derived optical thickness and LWP 
shows that the slope is 9.57 gm -2 and the intercept is 
— 1.35 gm~ 2 with the correlation coefficient of 0.88. If 
we omit the three large LWP values (>20gm~ 2 ) in the 
regression analysis, then the slope and the intercept of 
the regression are 5.0 and 0.16 gm~ 2 , respectively, and 
the correlation coefficient is the same (0.88) as that 
obtained including the three large LWP values. We 
mention that the peak values of the optical thickness 
derived with the lidar were consistently lower than the 
AERI ones. The reason for this is not likely to be the 
multiple scattering effects on the Raman nitrogen signal Jyj 
that can result in underestimation of the optical thick- 
ness because we obtained almost the same values from — I 
( 3 „ assuming the lidar ratio of 20 sr. C 

Figure 10 shows a scatterplot of LWP obtained with ^ 

the Raman lidar and AERI on 1 May; 5 June; 13 and ^ 

26 July; and 9, 10, and 31 August 2006 — on all of these 
days the overhead cloud was determined to be warm and 
could be fully profiled by the Raman lidar. The method ^ 
for computing LWP is the same as before. We selected 
the data for the cloud optical thickness smaller than 3 to 
ensure that the lidar fully profiled the cloud. The num- 
ber of data points used was 147. The correlation co- 
efficient between the data was 0.66. One probable 
reason for the low correlation is that the two instruments 
measured different cloud volumes because of 1) the 
difference in the location (—350 m) between the in- 
struments, 2) the difference in the FOV (0.2 mrad and 
the beam divergence of <0.1 mrad for the lidar and 
46 mrad for AERI, corresponding to the horizontal di- 
ameter of <43 cm for the lidar and 138 m for AERI at an 


altitude of 3 km), and 3) the inhomogeneity of the cloud 
distribution. We estimated the inhomogeneity of the 
cloud distribution from the high-temporal-resolution 
AERI retrievals (19-25-s resolution) and found that the 
standard deviation of the LWP were 10.6 and 4.1 gm“ 2 
for the period of Figs. 10 and 6, respectively. This sug- 
gests a larger cloud inhomogeneity in the data shown in 
Fig. 10 than those in Fig. 6 and helps to explain the better 
agreement shown in Fig. 6. Thus, for the validation of 



0 12 3 

Opt.Thick._Raman Lidar 

Fig. 9. Scatterplot of LWP as a function of optical thickness 
obtained with the Raman lidar for the same time period as in Figs. 5 
and 8. The equation in the figure is the result of linear regression 
with y = LWP, x = optical thickness, and r is correlation coefficient. 
Close-up for the small optical thickness values is the inset. 
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Fig. 10. Scatterplot of LWP for the values with the Raman lidar 
and AERI for cloud optical thickness lower than 3 for 1 May- 
31 Aug 2006 over ARM SGP site. Error bars show the standard 
deviation for AERI-derived LWP and measurement uncertainty 
estimated from the uncertainty in the lidar signal for the lidar- 
derived LWP. The equation in the figure is the result of linear re- 
gression with y = LWP(Raman lidar), y = LWP(AERI), and r is 
the correlation coefficient. 

the lidar measurements, it is better to compare mea- 
surements under homogeneous cloud conditions to min- 
imize the sampling differences between the instruments. 
We will show examples of such comparisons in the next 
subsection. 

c. Thick homogeneous clouds (low stratus) 
measurements over ARM SGP 

To minimize the influence of the sampling differences 
between the instruments, we analyzed homogeneous 
stratus clouds over the ARM SGP site in 2008. In this 
case, we compared the lidar-derived LWC with those 
obtained by MMCR-WVR, as we could not obtain 
LWP from the lidar data because the cloud was so 
thick that it strongly attenuated the laser beam and 
prevented full penetration of the cloud. The MMCR 
measures the vertical distributions of the reflectivity, 
vertical velocity, and Doppler spectrum width at 34.86 GHz 
(wavelength 8.66 mm) (Moran et al. 1998). The beam- 
width is 0.2° (3.5 mrad). The vertical distribution of 
LWC is derived from the MMCR-derived reflectivity 
and WVR-derived LWP using the method described 
by Dong et al. (1998) and Dong and Mace (2003). The 
MMCR-WVR-derived data were obtained from Cloud 
Retrieval Ensemble Dataset (CRED; Zhao et al. 2012) in 
the ARM Data Archive (http://www.arm.gov/data/evaP49). 


The resolution of the data was 1 h in time and 45 m in the 
vertical. The uncertainty in the MMCR- WVR-re trie ved 
LWC is 30% under the best conditions (Ebell et al. 
2010 ). 

Figure 11 shows the vertical distributions of LWC 
obtained with the lidar and MMCR-WVR for the pe- 
riod from 0300 to 0600 UTC 9 April 2008. The lidar data 
are averaged every 1 h (0300-0400, 0400-0500, and 
0500-0600). The derived vertical profiles of LWC from 
the lidar are consistently subadiabatic from the cloud 
base to 150 m above the cloud base (0.9-1.05 km in al- 
titude), demonstrating physical feasibility of the results, 
and agree generally well with the MMCR-WVR re- 
trievals. The calibration coefficients obtained by fitting 
the lidar data to the MMCR-WVR data in that altitude 
region were /lw = 0.12, Cwv = 0.03, and/ A E = 1.8, which 
are different from those obtained by comparison with 
AERI. The reason for the differences is possibly the 
measurement uncertainties of 30% in LWC for MMCR- 
WVR data (Ebell et al. 2010) and —10% for the lidar data 
estimated from the uncertainty in the lidar signals. The 
fractional contributions of the Raman water vapor to the 
LW channel and the aerosol fluorescence to the LWR 
are 30% and <1%, respectively, indicating that LWC is 
not sensitive to / AE for this case. 

Above 150 m from the cloud base, the laser could not 
penetrate effectively into the cloud but the WVR and 
MMCR radiowave of could penetrate it. The penetra- 
tion optical depth at which the lidar and MMCR-WVR 
results tended to diverge was approximately 3. This was 
estimated from /3 P assuming the lidar ratio of 20 sr. The 
total liquid cloud optical thicknesses obtained with 
MMCR were 11, 18, and 13 for those periods. 

We mention that the presence of signal-induced noise 
has been observed above 1.1 km in altitude and influences 
the results. Studies of the detectors’ response (not shown) 
indicate that detector nonlinearities may also be influ- 
encing the results shown for large optical depths. 

Figure 12 shows the temporal and vertical cross sec- 
tion of (3 p and LWC obtained from the lidar. The tem- 
poral resolution is 5 min and the vertical resolution is 
45 m. We can see in the figure the cloud base is stably 
present between 0.9 and 1.0 km for three hours. The 
value of LWC increased from 0 to 0.1 gm~ 3 from the 
cloud base to 150 m above that. 

Figure 13 shows the scatterplot of LWC obtained with 
the lidar and MMCR-WVR near the cloud base for the 
17 data points on 9 and 17 April 2008. Although the 
number of the data points is small, the correlation co- 
efficient between these data (0.82) is higher than that 
obtained with the comparison with AERI (0.66), partly 
because the homogeneity of the cloud distribution. How- 
ever, it would be preferable to increase the number of 
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Fig. 11. Vertical distributions of hourly averaged LWC obtained with the Raman lidar (black line with solid diamond) and MMCR- 
WVR (gray line with asterisk) for (from left to right) the period 0300-0600 UTC 9 Apr 2008 over ARM SGP site. Black dotted lines show 
the adiabatic profile. Error bars show the standard deviation for the MMCR-WVR and the measurement uncertainty estimated from the 
uncertainty in the lidar signals for the Raman lidar. 


comparisons in the future to increase the statistical confi- 
dence of the correlation. 


4. Discussion 

a. Aerosol fluorescence 

The results shown in the previous section (section 3c) 
were likely under the conditions of small influence of 
aerosol fluorescence on the liquid water signals because 
we did not observe any significant increase in the liquid 
water signal in cloud-free air. However, we have found 
some data that showed a large influence possibly due 
to aerosol fluorescence. Figure 14 shows the vertical 
distributions of LWC obtained with the Raman lidar 
and MMCR-WVR (left panel) and (3 p and depolari- 
zation ratio obtained with the lidar (right panel) for 
0800-0900 UTC 10 August 2006 over the ARM SGP 
site. It showed a peak in LWC with 0.025 gm~ 3 at an 
altitude of 3.9 km where /3 P showed a large peak with 


(a) Particle Backscattering Coeffient 



Time-UTC [Hour] 


0.055 knU 1 sr -1 . The MMCR-WVR LWC data also 
showed a peak at the altitude. This peak was probably 
cloud because the relative humidity computed from the 
lidar-derived water vapor mixing ratio and radiosonde - 
derived temperature and pressure was saturated with 
respect to water. The figure also shows a large peak in 
LWC with 0.4 gnC 3 at an altitude of 2 km where j3 p was 
0.0018 km - 1 sr -1 . This peak was probably not Raman 
backscattering from a water cloud because the relative 
humidity was unsaturated (55%). We suspect that this 
was the fluorescence by aerosols at the LW wavelength 
(403 nm). However, we have no information about the 
aerosol composition in that region for the period. The 
depolarization ratios are lower than 10% in that region, 
suggesting that spherical or near-spherical particles were 
predominant. The backward trajectory computed by 
Hybrid Single Particle Lagrangian Integrated Trajectory 
(HYSPLIT; Draxler and Rolph 2012) have shown that 
the air parcel of the aerosol layer had been carried over 
the southern part of United States for 3 days before, 


(b) Liquid Water Content 



Time-UTC [Hour] 


Fig. 12. Temporal and vertical cross sections of (a) (3 P and (b) LWC obtained with the Raman lidar for the period 0300-0630 UTC 9 Apr 

2008 over ARM SGP site. 
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Fig. 13. Scatterplot of LWC obtained with the Raman lidar and 
MMCR-WVR on 9 and 17 Apr 2008 over ARM SGP site. Error 
bars show the standard deviation for the MMCR-WVR and 
measurement uncertainty estimated from the uncertainty in the 
lidar signals for the Raman lidar. The equation in the figure is the 
result of linear regression with y = LWC(Raman lidar), x = LWC 
(MMCR-WVR), and r is the correlation coefficient. 

when the presence of smoke in Mississippi and Alabama 
was predicted by the global aerosol model (Navy 
Aerosol Analysis and Prediction System model, de- 
scription and results available from the web pages of 
the Naval Research Laboratory, Monterey, California; 


http://www.nrlmry.navy.mil/aerosol). This suggests that 
smoke particles were the source of the fluorescence in 
the aerosol layer. However, we cannot make this conclu- 
sion with certainty because many atmospheric aerosols 
exhibit fluorescence: chemically aged terpenes, fungal 
spores, polycyclic aromatic hydrocarbons, bacteria, cel- 
lulose pollens, and mixtures of various organic carbon 
compounds (Pinnick et al. 2013; Lighthart and Mohr 1994; 
Jaenicke et al. 2007). We found that the lidar data in 
August 2006 showed a similar increase in LWC in the 
lower and midtroposphere in the absence of clouds. To 
confirm that aerosol fluorescence can interfere with the 
Raman LW signal, we measured the fluorescence spec- 
trum of pine pollen in the laboratory at GSFC. To measure 
the fluorescence spectrum, we used a spectrofluorome- 
ter (Fluorolog, HOR1BA Jobin Yvon). The incident 
light wavelength was 354.7 nm with 1-nm width and the 
detection angle is 167.5° to the incident light. Figure 15 
shows the spectrum of the pollen fluorescence, which 
has a broad increase at wavelengths larger than 365 nm 
that overlaps the Raman liquid water scattering band and 
the transmission curve of the LW channel. 

Although the absolute fluorescence cross section was 
unavailable in these laboratory experiments, we can 
expect that the aerosol fluorescence may interfere with 
the LW signal under the conditions of high fluorescence 
aerosol concentration in the atmosphere. For example, 
Gelbwachs and Birnbaum (1973) reported the aerosol 
fluorescence intensity at the wavelength range of 610- 
660 nm excited at 488 nm in the ambient atmosphere 
on a smoggy day in California was 7 times as large as that 
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Fig. 14. (a) Vertical distribution of LWC obtained with the Raman lidar (black) and MMCR-WVR 
(gray) for 0800-0900 UTC 10 Aug 2006 over ARM SGP site, (b) Vertical distributions of particle 
backscattering coefficient (solid) and depolarization ratio (dotted) obtained with Raman lidar for the 
same time period. 
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Fig. 15. Fluorescence spectrum of pine pollen (thick solid line.) 
Thin solid line shows the Raman spectrum of liquid water. Thin 
dotted line shows transmission of IF for the LW channel of the 
ARM lidar. 


due to N0 2 at 0.25 ppm. This intensity (2.4 X 10 -11 m _1 
sr” 1 ) would correspond to 0.2 g m 1 in LWC using the 
parameters given by Gelbwachs (1973) of the N0 2 ab- 
sorption coefficient of 10 -3 m -1 ppm -1 at 488 nm, the 
factor by which the N0 2 fluorescence quenched at 1 atm 
of 2 X 10~ 5 and detection wavelength width of 6.5 nm 
for the ARM lidar (Table 2). This result indicates that 
fluorescence is very likely to be a contributor to the 
liquid water signal under the conditions observed by 
Gelbwachs and Birnbaum (1973). Immler et al. (2005) 
have also observed inelastic scattering from an aerosol 
layer in the lowermost stratosphere with the Raman 
lidar at a wavelength at 407 nm with a bandwidth of 
5.6 nm. They attributed it to the fluorescence by smoke 
particles from a forest fire. Sugimoto et al. (2012) have 
observed broad fluorescence at wavelengths shorter 
than 510 nm excited at 354.7 nm from Asian dust and air- 
pollution aerosols transported from urban and industrial 
areas in the lower troposphere. The influence of fluo- 
rescence on the WV signal is generally less because of 
the typical use of narrow bandwidth filters (e.g., 0.3 nm) 
for water vapor detection. However, care must be taken 
in analyzing Raman water vapor lidar data when mea- 
suring very low water vapor concentrations under condi- 
tions where aerosols may be present because fluorescence 
may be artificially enhancing the signal even for lidar 
systems that use narrow band detection. 

b Liquid water Raman backscattering cross section 

We estimate the Raman backscattering cross section 
of liquid water from / LW that was determined by the 
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Fig. 16. Raman backscattering cross section of LW at an exciting 
wavelength of 354.7 nm estimated from this study (solid diamonds) 
and those reported by laboratory experiments (asterisk). Error 
bars of this study indicate the uncertainties estimated from the 
measurement uncertainties of the two instruments. Error bar of 
Faris and Copeland (1997) indicates the single standard deviation 
of the measurement. Values of Li and Meyers (1990) and Ahmad 
and lies (2001) are overlapped. 


comparisons of the lidar data with AERI and MMCR- 
WVR (sections 3b and 3c). In addition, we show the 
result estimated using the data obtained between 1.1 
and 1.5 km for 0800-0859 UTC 3 May 2008, when the 
correlation coefficient between the lidar-derived LWC 
and MMCR-WVR-derived LWC is the highest (r = 
0.96) with/ L w = 0.18, Cwv = 0.02, and/ A E — 0.4. To 
estimate the liquid water Raman cross section from /| W , 
we used Eq. (2) and the parameters given in section 2d 
(dawvIdLl, (wv. £lw, and G) except for the liquid water 
cross section. Figure 16 shows the results as well as those 
reported from recent laboratory experiments. The uncer- 
tainties are calculated from the measurement uncertainties 
in LWC of each instrument and the uncertainty in the 
water vapor Raman cross section using standard error 
propagation. They are 4% for AERI (Turner 2007), 
30% for MMCR-WVR (Ebell et al. 2010), 18% in 
August 2006 and 5% in April-May 2008 for the Raman 
lidar estimated from the uncertainty in the lidar signals, 
5% in C W v (section 2d), and 10% in the water vapor 
Raman cross section based on Penny and Lapp (1976). 
The Raman liquid water cross-sectional values estimated 
from the comparisons vary from (2.3 ± 0.7) X 10 -33 to 
(3.4 ± 1.1) X 10 _33 m 2 sr _1 with an ensemble mean of 
2.7 X 10 -33 m 2 sr -1 , which is 46% lower than those of 
Ahmad and lies (2001) and Li and Meyers (1990), and 
28% lower than that of Faris and Copeland (1997). The 
relative standard deviation of the ensemble mean for the 
lidar data is 23%, which is close to the maximum un- 
certainty of the cross section reported by Faris and 
Copeland (1997). A reasonable explanation for the 
discrepancy is that errors may have been introduced in 
the detection efficiency ratio of LW and WV channels, 
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which we assumed to be one. As mentioned earlier, the 
sensitivity of individual PMTs can differ significantly. It 
would be possible to measure the ratio of detection ef- 
ficiency of these two channels by interchanging the 
PMTs or IFs (Vaughan et al. 1988). The lamp mapping 
technique is also a promising solution for determining 
the overall detection efficiency of the receiving system 
(Venable et al. 2011). However, many more lidar cases 
are needed to increase the number of comparisons to 
increase the statistical confidence of the estimate of 
Raman cross-section value because the large differences 
of /'i vv- C wv , and/ AE in thin and thick clouds (Table 3) 
could be due to the differences in the cloud character- 
istics (i.e., LWC, droplet size distribution, and number 
concentration) and the aerosol fluorescence effect. 

5. Summary and conclusions 

This paper has described the current state of our re- 
search into the Raman liquid water cloud measure- 
ment technique based on the results of spectroscopic 
measurements acquired at NASA GSFC and the anal- 
ysis of data from the ARM SGP site. The conclusions of 
this study are as follows: 

1) The observed backscattering spectra in/below cloud 
were spectrally similar to laboratory and theoretical 
calculations of Raman liquid water and water vapor 
spectra over GSFC on 15 November 2011. The 
intensity ratio of liquid water to water vapor in the 
cloud was much higher than that below the cloud and 
was consistent with reasonable values of cloud liquid 
water content. 

2) The coefficients for computing LWR from the ARM 
lidar were determined by comparison with AERI 
and MMCR-WVR to be / LW = 0.12-0.18, C W v = 
0.01-0.03, and/ AE = 0.4-1. 8 for the data in 2006 and 
2008. The possible reason for the large ranges in 
these coefficients is the measurement errors of the 
instruments (4% for AERI, 30% for MMCR-WVR, 
and 2%-18% for the Raman lidar) and the sampling 
differences of the two instruments. 

3) Strong qualitative correlation was obtained in the 
LWP time series of lidar and AERI for some cases, 
although considerable scatter in the overall ensem- 
ble of comparisons may be due to the different fields 
of view of the instruments, distance between the 
instruments, and the horizontal inhomogeneity of 
the clouds. 

4) The LWC of a homogeneous low-level stratus cloud 
was obtained near the cloud base to 150 m above 
where the optical depth was less than 3 and was in 
good agreement with retrievals from MMCR-WVR. 


5) We show a significant enhancement of the liquid 
water signal in the lidar measurement in August in 
2006 that we believe to be due to aerosol fluores- 
cence. This is a concern both for liquid water and 
water vapor measurements by lidar. 

6) The Raman backscattering cross section of liquid 
water obtained from this study is 28 %— 46% lower 
than the most current laboratory values for the 
Raman liquid water cross section, although possible 
differences in lidar PMT detector efficiencies have 
not yet been accounted for. The uncertainty of the 
Raman lidar results nearly overlaps with those of 
Faris and Copeland (1997). 

Suggestions for improving the Raman lidar liquid 
water cloud measurement technique are as follows: 

1) Measure the cloud spectrum with high S/N ratio 
under the conditions of high/low aerosol loading to 
study the influence of the aerosol fluorescence on the 
Raman liquid water channel. 

2) Increase the number of comparisons of the lidar- 
derived cloud properties with those obtained with 
the other instruments (e.g., AERI, MMCR-WVR, 
and aircraft-based instruments) to further test the 
lidar measurement. 

3) Measure the lidar detector’s efficiency ratio of LW 
and WV (or N 2 ) channels to reduce the uncertainty 
in the Raman backscattering cross section of liquid 
water that is estimated from the comparison. 

4) Work on correcting for the nonlinear response of the 
detectors in the presence of thick clouds and/or to 
improve the linearity of the detector’s response to 
the strong signal. 

With success at reducing the measurement uncer- 
tainties of the Raman lidar cloud liquid water retrievals, 
the Raman lidar can potentially provide a complete set of 
measurements for studying the aerosol indirect effect 
without concern for the different fields of view, and tem- 
poral and spatial resolutions of the different instruments. 
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